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Abstract 

We analyze the charge-parity (CP) asymmetry in the charged Higgs boson decays to chargino- 
neutralino pairs, — > xfx^^ i = 152, j = 1,...4. We show first that these modes have a 
large branching ratio for mfj± > 600 GeV. We use Cutkosky rules to obtain the analytical for- 
mulas needed for the evaluation of the asymmetry under consideration. We then calculate the 
CP asymmetry in chargino-neutralino decays by including supersymmetric mass bounds, as well 
as constraints from h — > 57, {g — 2)^, Ap and electric dipole moments. Finally, we discuss ob- 
servability of the asymmetry at the LHC by calculating the number of required charged Higgs 
events to observe the asymmetry for each decay channel. We show that the inclusion of constraints 
considerably reduces the projected CP asymmetry, and that the optimal channel for observing the 
asymmetry is X2- 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) at CERN is expected to unravel the mysteries of 
electroweak symmetry breaking as well as provide signals of physics beyond the Standard 
Model (SM). One would hope that in the first instance, the Higgs boson would be found. In 
the SM, there is only one physical (neutral) Higgs boson. Should a charged scalar particle 
be found, one would have to unravel the underlying physics, as most frameworks beyond SM 
predict charged Higgs bosons. The simplest and the most popular of such models are the 
two Higgs doublet model(s) and supersymmetric models. One could distinguish among the 
"new Higgs" by testing their charge-parity (CP) properties. In the (SM) the CP breaking 
occurs through only one weak phase in Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 



observed direct CP violations in K 



l| and B 



^ decays can be accommodated via the CKM 



matrix in the SM. Models beyond the SM contain several new sources of CP violation, some 
of which could affect Higgs boson decays. The minimal supersymmetric standard model 
(MSSM) contains many parameters which can in principle be complex, even after making 
all allowed rotations to get rid of unphysical phases. 

Such new sources of flavor and CP violation give rise to enhanced CP violation, some 
of which could provide distinguishing signs for MSSM at present and future colliders. In 
a previous work [sl, we discussed the CP asymmetry in charged Higgs decays Uidj 
in the framework of the MSSM, where Uidj = cb,cs,ib,is. We showed that, although the 
channels is and cs have sizable CP asymmetry with respect to the channels ib and cb, this 
result has to be taken with caution since the formers have very small branching ratios (i?r's) 



which makes them harder to observe. We calculated {A'^p x Br)~^ {4] for each decay mode, 
which is proportional to the number of required charged Higgs bosons produced at colliders^ 
to observe an asymmetry for a given channel 0, [gJ . Based on this analysis, we concluded 
that H~ cb and H~ — > tb are the the optimal channels which could reveal a measurable 
CP asymmetry at the order of 10 — 15%. 

In this study, we calculate the CP asymmetry in charged Higgs decays in the purely 
supersymmetric mode H"^ xfxp = 1)2, j = 1,...4, compare the size of the CP 



^ {A'^,p X Br)^^ is closely related to the total number of events N required to establish a measurable CP 
violation for a particular mode. The exact formula is = s^(A^p x Bre)^^. Here s is the standard 
deviation and e is the detection efficiency. 
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asymmetry for each i and j , and evaluate the number of charged Higgs needed to observe each 
asymmetry. We show that the asymmetry can be very large if we include only constraints 
requiring supersymmetric particles to have masses consistent with experimental bounds. 
However, if we also include low energy phenomenology constraints coming from b 57, {g — 
2)^, Ap and electric dipole moments, the predicted asymmetry and the allowed parameter 
space are considerably reduced. 

The outline of the paper is as follows. In the next section, we briefly review the CP 
structure of the MSSM. The processes and the method of calculation are outlined in Section 
mil and the numerical analysis of the decays under consideration are presented in the Section 
HVl Finally, we conclude in Section |Vl and give the relevant analytical expressions in the 
Appendix. 



II. THE CP STRUCTURE IN THE CHARGING AND NEUTRALINO SECTORS 
OF THE MSSM 

In addition to the CP violation induced by the CKM matrix in the quark charged current, 
CP phases can be introduced explicitly in the MSSM, by complex Yukawa couplings of Higgs 
bosons to quarks and squarks. New supersymmetric-only sources of CP violation in MSSM 
are: the soft supersymmetry breaking gaugino masses M^, a = 1,2, 3; the bilinear Higgsino 
coupling parameter /i; and the soft supersymmetry breaking trilinear scalar coupling of the 
Higgs bosons with scalar fermions f , Af, f = u, d, c, s, t, b, e, /i, r. In principle, each of these 
parameters can have independent CP-phases, making a general analysis of CP violation or 
asymmetry very complicated. We make several simplifying assumptions. First, we assume 
that at the unification scale the gaugino masses have a common phase and the trilinear 
couplings are all equal and have another common phase. In order to avoid known problems 
with the Electric Dipole Moments (EDMs) one could deviate from exact universality 
and consider A f to he diagonal in flavour space with vanishing first and second generation 
couplings. This leaves independent phases in /i, and Af. However, the symmetries of the 
MSSM can be used to re-phase one of the Higgs doublet fields and the gaugino fields such 



that Ma are real 



In addition, the CP-violating phases associated with the sfermions 
of the first and second generations are severely constrained by bounds on the EDMs of the 
electron, neutron and muon 10| and this limits arg(/i), which cannot exceed values of the 



3 



order of lO'^ - IQ-^. 

There have been several proposals 



ll|-[l3| to evade these constraints without suppressing 



the CP-violating phases. One possibility is to arrange for partial cancellations among various 



contributions to the EDMs [13[. Another option is to make the first two generations of 
scalar fermions rather heavy, of order 3 TeV, so that the one-loop EDM constraints are 
automatically avoided. One could invoke the effective SUSY models [3] where decoupling 
of the first and second generation sfermions are used to solve the SUSY Flavour Changing 
Neutral Current (FCNC) and CP problems without spoiling the naturalness condition. We 
adopt the latter version of a CP- violating MSSM for our analysis, along with Af = for the 
first two generation sfermions. We also simply neglect arg(/i) and consider fi real. Thus, 
we assume that the only non-zero phases emerge from the trilinear couplings At h T and we 
assume a common phase^ aig{At) = arg(ylf,) = aig{Ar) = arg(yl) unless otherwise stated. 

The superpotential W of the MSSM Lagrangian and the relevant part of the soft breaking 
Lagrangian L^^ft respectively 

W = fiH^H^ + Yi'H'Ve^ + Y^^'H^Q'd^ + Y^H^Q'u'j^ (2.1) 

where H^, and are the Higgs doublets with vacuum expectation values vi and V2 re- 
spectively, Q is the SU{2) scalar doublet, m, d are the up- and down-quark SU{2) singlets, 
respectively, Q,u,d represent scalar quarks, Yu^^ are the Yukawa couplings and i,j are gen- 
eration indices. Here A^^ represent the trilinear scalar couplings. 

The CP violation effects enter the Higgs decays to chargino and neutralino at one loop 
only through loops containing scalar quarks. For simplicity, and to avoid introducing new 
parameters, we neglect flavor mixing in the scalar quark mass matrix. The scalar quark 
mass is taken as: 

■Wk= I "" 1 (2-3) 



The CP-violating phases arg(/i) and arg(Af_f,) could in principle be measured directly in the production 



cross sections and decay widths of (sVparticles in high energy colliders [8|, IJ]- 16[ or indirectly via their 



radiative effects on the Higgs sector [8, 15|. 
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with 



= M|_^ + m2 + cos2/3(T,-Q,4)M|, 
At = At - fi cot (3 , Ab = Ab- fitanf] . (2.4) 

The mass squared matrix for the scalar r lepton is similar to the one for b-squark, with the 
appropriate substitutions. Here we assume fi real and a common phase for At and Ab. Note 
that we neglect the effects of Ac^g since they are multiplied by the charm or strange quark 
mass. Here tan/5 is the ratio of the vacuum expectation values of the two neutral Higgs 
bosons. 

We summarize the basic properties of the charginos and neutralinos entering in the cal- 
culation of the CP asymmetries, as well as give the basic couplings between them and the 



scalar quarks entering the loop calculations. For more information, see [17|. 

The charginos xt ~ 2) are four component Dirac fermions. The mass eigenstates 
are mixtures of the winos and the charged higgsinos H^, H2 ■ The chargino masses are 
eigenvalues of the diagonalized mass matrix, X, given by 

/ M2 \/2Mvi/sin/5 

X = 

\ ^/2Mw cos/5 /i 

where M2 is the soft SUSY-breaking SU{2)l gaugino mass parameter and /i is the Hig- 
gsino mass parameter. The relevant Lagrangian terms for the quark-down squark- chargino 
interaction are given by 

2 3 



(2.5) 



^nJx+=E E (Yt' KckmU Pl xt {^dY^ I - [g U,, [KckmU Pr xt {^dT da 

a=l i,j=l 

+ [U„2 {KcKM Yt%,] Pr xt d)'^'^''^'' 4} + H.c. , (2.6) 



where the index a refers to chargino mass eigenstates. Y^^2^ are the diagonal up- and down- 
quark Yukawa couplings, and V ^ U the chargino rotation matrices defined by U*^V~^ = 
diag(m^+,m^+). 

Neutralinos Xi (^ = 1;2,3,4) are four-component Majorana fermions. The mass eigen- 
states are mixtures of the photino, 7, the zino, Z, and the neutral higgsinos. Hi and if^ ■ 
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Their mass matrix Y is given by 

^ Ml -Mz sin cos /3 sin ^vi/ sin /3 ^ 

M2 Mz cos 6w cos [3 —Mz cos 6w sin [3 

—MzsmOw cos (3 Mzcos6wcos(3 — /x 

^ Mz sin 6w sin /5 — cos 6w sin /? — yU ^ 

The relevant Lagrangian terms for the quark-up squark neutrahno interaction are 

^ 9 n n _ .0 ^-n \(i+3)a ~ - at* i^diae n _ .0 ^"n Nia^, 

(*+3)a~ 



(2.7) 



n=l 1=1 ^ * 

- (^Nn2 + ^iVni tan Ow^ Pr X° (Ff;)*'^ - 7V„4 1;'''"^ Pfl Xn i'Tu) 



(2.8) 



where N is the rotation matrix which diagonahzes the neutrahno mass matrix M^o, 
N*Y N^^ = diagfm^o, m^o, m^o, m,^o). 

Xl ' X2 ' X3 ' X4' 

In aU of the above interactions, CP violation arise from aig{At h). Though we include 
squark contributions only for simplification, sleptons contribute too in an analogous manner. 
We check the sensitivity of the CP asymmetry to arg(yl^) as well. 

III. CHARGED HIGGS DECAYS xfx] AND CP ASYMMETRY 

The relevant Feynman diagrams for calculating the CP asymmetry in charged Higgs 
decays to chargino-neutralino pairs are shown in Fig. [H There are indeed more than 120 
Feynman diagrams contributing at one-loop level but based on our choice of the CP phases 
the list goes down to eleven diagrams given in Fig. [TJ 

The CP asymmetry is defined as: 



T{H-^X-x'j)+nH^^xtxl 



For simplicity, in the denominator of Eq. (13.11) . we approximate the decays width of the 
charged Higgs by their tree-level contribution since only the real part of the decay width is 
needed, and including one-loop contributions is not essential. 

The CP-odd observable Acp requires a nontrivial phase from Feynman diagrams (called 
absorptive or strong phase), in addition to the weak phase from arg(y4tb,-) mentioned 
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in the previous section. This way, the imaginary part of the amphtude is non-zero, 
Xm(Amphtude) 7^ 0. One way of introducing such a phase is through one-loop Feynman 
diagrams, where some of the intermediate particles go on-shell. As tree-level contributions 
do not generate a CP asymmetry by themselves, we estimate the CP asymmetry from the 
interference of tree and one-loop diagrams, which dominate over corresponding contributions 
from loop-loop terms. We show in Fig. [1] the tree level and the one-loop graphs which only 
contribute to the asymmetry [i.e., the ones with scalar fermions in the loop). 

One convenient way of calculating CP asymmetry is to extract the absorptive part of the 
loop diagrams by applying the Cutkosky rules 13]. For self-energy type diagrams there is 
one possible cut, but for the vertex diagrams there are three types of cuts depending on the 
masses of the external particles. So, in Fig. [21 we show possible cuts where the internal-cut 
states represent all possibilities in each case. For the details of the method of calculation. 
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FIG. 1: The tree and relevant one- loop diagrams contributing CP asymmetry in t' Hooft-Feynman 
gauge for the decays —>■ XcXn- 
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FIG. 2: The unitarity cuts for self energy (a) and vertex type (6) generic diagrams. 



see Ref. 

For self-energy diagrams the vertical one can be: a cut through scalar top - down quark, 
scalar bottom - up quark, or scalar top - scalar bottom loop. They contribute if m^- ^ 

+ ^dji fTT'^- ^ + f^uj, or rrifj- ^ + m^, respectively. Here Xc represents the final 
state chargino with c = 1,2. 

For the vertex diagrams, as seen from Fig. [21 there are three possible cuts: one vertical, 
two horizontal ones. The vertical cuts through r — z/^, b — t, and t — b always contribute if 
a non-zero phase for A^, At, and Ai, is assumed, respectively, and for charged Higgs masses 
> rrit. The other vertical cuts through f — Ur and t — b depend on the chosen Higgs 
mass as well as the phases of At^b,T- The situation is quite similar for the other two horizontal 
cuts, but this time the mass of the chargino or neutralino has to be greater than the sum 
of the masses of the internal particles where we cut through. Some analytical formulas are 
given in Appendix for the horizontal and vertical cuts in a generic way, rather than giving 
lengthy analytical results for the imaginary parts. 

We generate the amplitudes using FeynArts and proceed with our own code for the 
rest of the calculation. In the numerical calculation we checked our results for the imaginary 
parts of the two-point and three point diagrams with the package program LoopTools 19 1 
and found close agreement. Also, for the numerical diagonalizations of chargino, neutralino 
and scalar quark mass matrices, we cross checked our routines with the Hahn's routines 
Diag Q. 
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IV. NUMERICAL ANALYSIS 



In this section we present a comparative numerical analysis of the CP asymmetry in 
charged Higgs decays into all possible pairs of chargino and neutralino states. Our aim is 
two-fold: the first is to predict the CP asymmetry and indicate the optimal channel as well 
as estimate the number of charged Higgs needed to probe it. The second is to analyze the 
effect of various experimental constraints such as b ^ 57, {g — 2)^, Ap and EDMs and 
determine their significance. 

As we approximate the partial decay widths of the channels by their tree level contribution 
in the denominator of Eq. (13.11) . the denominator becomes 2r'^(if~ X7X^) where r° 
represents the tree level partial width. The partial decay widths of the decay into all possible 
eight channels of chargino-neutralino pairs with tree level approximation are shown in Fig. [3] 
as a function of the Higgs mass tuh-, for some values of the free parameters of MSSM (the 
common scalar mass scale M^vsy, the gaugino mass parameter M2, the trilinear coupling Af, 
fi, and tan j3). The decay widths are plotted on the left for intermediate tan (3 value, while the 
ones on the right are for large tan (3. For the calculation of the rest of charged Higgs channels, 
the FeynHiggs program 2l| is used. Among these, there are fermion channels; H~ fiUf 
among which rur is the main decay mode for low Higgs mass region {rriH- < 140 GeV) 
and tb is for intermediate Higgs mass values {rriH- > rrit), the Higgs- vector boson channels 
— s> h^W, H^W, A^W, and the sfermion channels — ^ fifj i,j = 1...3 where they are 
usually significant for heavy Higgs mass region. 

As seen from Fig. [31 at intermediate tan/9, some of chargino-neutralino channels start 

being competitive and dominant with respect to the sum of the other open channels for a 

Higgs mass around rriH- ^ 600 GeV. This is not the case at large tan/5 (the case for low 

tan/3 is similar to the case for high tan/3 and we do not graph it separately). In these cases, 

in the heavy Higgs mass region, the eight chargino-neutralino channels split into two groups; 

H~ XiX2^ X-iXi^ H~ X-iXl with the branching ratios less than 0.006 and 

the other five channels with the Br's around 0.1-0.3. Since usually a smaller branching ratio 

fl 

is an indication of larger observable CP asymmetries (6|], we could expect that H ^ Xi X25 
H~ — s> X2X1J cind — s> X2XA may yield observable CP asymmetries. Of course an extensive 
numerical analysis is required and in the rest of the section we investigate this by scanning 
the parameter space. 
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FIG. 3: On the left: MsusY = 500 GeV, fi = 250 GeV, tan/3 = 10, Ma = 200 GeV, and \ At^b,T\ 
400 GeV. On the right: Same as on the left but for tan /3 = 30. 



We first check the sensitivity of the CP asymmetry to the free parameters and phases 
and then scan the sensitive parameters for the maximal asymmetry on a multidimensional 
parameter space. At each stage we show both the regions allowed by imposing only su- 
persymmetric mass limits, and then show the parameter regions excluded by imposing the 
experimental constraints coming from b 57, {g — 2)^, Ap and EDMs. The current exper- 
imental mass lower bounds considered are mu,f > 96 GeV, > 89 GeV, m^o > 50 GeV, 



and > 94 GeV 
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FIG. 4: On the left: MsusY = 1 TeV, /x = 250 GeV, tan/3 = 10, M2 = 200 GeV, \At,b,T\ = 1-2 
TeV, aTg{Ah) = arg{AT-) = arg(Af) = arg(A) and m^j- = 1 TeV. On the right: Same as on the left 
but arg(A) = tt/2. Excluded regions are disfavored by all experimental bounds considered here. 



The constraints we included to restrict the CP asymmetry are: {g — 2)^j < 40 x 10 23 | 



and |Ap| < 10 ^ 22|; as well as electric dipole moment constraints: \d^\ < 6 x 10 e.cm 



for neutron 
mercury 



2^, I^Thl < 1-6 X 10- 



-27 



e.cm 



25|. For 6 — i> 57 we use at 3a [26| 



or thallium^ 



2J], and |ciHg| < 2 x 10"^^ e.cm for 



2.53 X 10"^ < Br{b 57) < 4.34 x 10" 



(4.1) 



The calculation and inclusion of these bounds have been done with the use of FeynHiggs 



2l|. We don't impose SUSY GUT relations for gaugino masses. Since the 6 ^ 57 decay 
can be accommodated by a rather heavy gluino mass, and this does not influence our results 
for the CP asymmetry, this assumption allows us to keep M2 rather light. So, we assume a 
gluino mass of nig = 5 TeV but express Mi in terms of M2 using GUT relations. We find 
that inclusion of the experimental constraints not only restricts the parameter space but 
also has a significant effect on the CP asymmetry overall. 

In Fig. m the CP asymmetry is depicted as a function of a common trilinear coupling A 



The upper bound is indeed the translated one on the electron EDM de 
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FIG. 5: On the left: MsusY = 1 TeV, tan/? = 10, M2 = 200 GeV, \At^b,T\ = 1-2 TeV, 
arg{Ah) = arg(Ar) = aig^At) = tt/2 and mjj- = 1 TeV. Excluded region is disfavored by mass 
lower bounds. On the right: Same as on the left but = 200 GeV. Excluded region is disfavored 
by all experimental bounds considered here. 

for the third generation lepton and quarks, assuming they share a common phase, arg(y4). 
The Higgs mass is assumed rather heavy, 1 TeV. The maximal CP asymmetry is obtained at 
arg(A) = 7r/2 which is still allowed after excluding the experimentally disfavored region. The 
largest asymmetry is obtained for the channels H~ X2X\i H~ — > X1X21 ^^'^ H~ — > X2X\- 
As expected the larger |y4| is, the bigger the asymmetry becomes. In these figures, almost 
all the contribution to the asymmetry comes from the interference of the vertex diagram 
7^ 7 in Fig. [T]with tree level. The asymmetry at maximum is around 0.1% for H" X2XA 
and much smaller for the rest. Although we did not include it in a separate figure, we also 
analyzed the CP asymmetry as a function of the phase of Ab by assuming the other phases 
zero, and separately for the phase of A^- and found no significant effect. Assuming no phases 
from Ab and Ar leads to no contributions to Acp from the diagrams #3,6,8,9, and 10. 
This has been verified numerically. 

In the next figure. Fig. [5l as we vary fi in the (—500, 500) GeV range, the magnitude of 
the CP asymmetry becomes as large as (0.3 — 0.5)% for H~ — ^ XiX2y ^ ~^ X2X% cind 
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FIG. 6: For the first row: MsusY = 1 TeV, tan/3 = 10, |At,f,,^| = 1.2 TeV, arg(A) = arg(^^) = 
arg(^t) = 7r/2. Left: nifj- = 1 TeV. Right: M2 = 200 GeV. For the second row: Same as above 
but |j4j = 1 TeV. Excluded regions are disfavored by all experimental bounds considered here. 
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H~ — s> X2 Xs- This time the region exduded violates the experimental lower bounds on the 
lightest chargino and neutralino mentioned above^. The asymmetry is sensitive with respect 
to variations in tan/3 and reaches (3 — 4)% in the allowed region for (again) H~ — XiX2- 

The next figure, Fig. El presents the asymmetry as a function of M2,mH-, and Msusy- 
The asymmetry vanishes for all channels except H~ XiXi ^^'^ H~ — XiX^ larger M2 
values for kinematical reasons (the heavier channels cannot be produced in the decay). As 
a function of the Higgs mass, the CP asjTiimetries for XixH ^^"^ H~ — > X2X\ are the 

most sensitive. The EDM constraints especially allow neither Msusy nor nifj- to be light. 
Choosing a smaller A value with a smaller phase would relax this restriction significantly, 
but would also decrease the CP asymmetry. 

The final part of the numerical analysis is devoted to the scan for maximal CP asymmetry 
among the four of eight channels (two light from Fig. [7] and two heavy ones from Fig. [8]) 
over sensitive variables of the parameter space (Msusy, m//-, M2, /i, tan /3). We choose the 
following four channels; H" XiX2^ X2X1, H' X2 X3, and X2XI, which 

seem most promising to yield larger asymmetries. For the sake of efficiency, we included 
some of the results of the previous parameter dependence. Since we verified that the effect 
of the phases aig{Ah) and arg(yl^) are negligible, we set them to zero, while keeping the 
absolute value of Ah, since there are diagrams with scalar- up/scalar-down quark loops which 
are not negligible. We also fix the phase of At as 7r/2, as we know that it maximizes the 
CP asymmetry. The absolute values of At and At, are chosen as 1.2 TeV, which appear 
to be values in the region which maximizes the asymmetry. Then we vary the rest of the 
parameters in the following ranges; Msusy G (600, 1200) GeV, ttih- G (600, 1500) GeV, 
M2 G (200, 1000) GeV, /i E (-500,500) GeV and tan/3 G (2,50). 

In Figs. [7]and[8|, we present the scatter plots of the modes H~ XiX^ H~ — > X2X1) 
H" — > X2 Xs! and H~ — ^ ^2^X4 i^i the {A^p x Br)~^ —Acp plane. As before |3|, we introduce 
here the quantity {A^p x Br)~^ as a measure of the number of required charged Higgs bosons 
for the observibility of the asymmetry. This is an order of magnitude estimation and the 
exact value can be obtained by multiplying it with the factor s^/e where s is the standard 
deviation and e is the detection efficiency. There are two types of points on the scan, color- 

^ The boundary points of the excluded region vary with respect to each decay mode. In Fig.[Sl we consider 
some average values for all of them. 
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FIG. 7: The scatter plot for the charged Higgs decays H Xi(2)X2{i) ™ {A'^p x Br) ^ — Aqp 
plane. These are obtained by scanning the sensitive parameters A £ (0, 1400) GeV, mjj- £ 
(600, 1500) GeV, tan /3 G (1,50), Afs G (200, 1000) GeV, /i G (-500, 500) GeV, and MsusY e 
(600, 1200) GeV. The phase is fixed, arg(A() = 7r/2. The blue points satisfy only the mass bounds 
and the green ones satisfy all constraints. The x-axes are in the logarithmic scale. The small 
graphs inside each graph represent the positive asymmetry Acp in the logarithmic scale. 

coded as blue and green. The blue dots only obey the basic mass lower bounds. As one can 
see, under these requirements only, each decay mode can generate a CP asymmetry as large 
as 30%^. Of course the input set for the parameters is different for each channel. These are 
listed for the maximal asymmetry points in Table [B The number of charged Higgs bosons 
needed to observe the predicted CP asymmetries is in the range of (10^ — 10^) x s^/e. In 
the small panels in Figs. [7| and [H] we also showed the distribution of points in logarithmic 
Acp scale. 

However, if in addition to the mass bounds, one adds the experimental constraints on 
b S7, {g—2)fj_, Ap and EDMs, only the green points on the graphs survive and the maximal 



^ The mode H X2X1 can also have blue points as large as 30% but we restricted the extent of the 
vertical axis for better resolution. 
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FIG. 8: The scatter plot for the charged Higgs decays H X2 ^3(4) i^i the {A^p x Br) ^ — Acp 
plane. These are obtained by scanning the sensitive parameters A G (0, 1400) GeV, niu- £ 
(600, 1500) GeV, tan /? E (1,50), M2 G (200, 1000) GeV, /i G (-500, 500) GeV, and MsusY G 
(600, 1200) GeV. The phase is fixed, arg(^j) = 7r/2. The blue points satisfy only the mass bounds 
and the green ones satisfy all constraints. The x-axes are in the logarithmic scale. The small 
graphs inside each graph represent the positive asymmetry Aqp in the logarithmic scale. Thus, in 
the small panel for X2 X3) Acp is multiplied by -1. 

asymmetries for each channel go down to few percent level, with the requirement of a larger 
number of Higgs to be produced for better statistics. Again, the maximal points and the 
corresponding input set and the required Higgs bosons are given in Table [T] as separate 
columns for comparison. As seen from the Table, for maximal asymmetry all four channels 
favor a heavy charged Higgs mass greater than 1 TeV and also either large or small tan j3. In 
the case where we include only the mass bounds, the low MgusY values in the chosen range 
(600, 1200) GeV are favored since in such cases rriH- > niu. + m^, is satisfied, the squarks 
can go on-shell in the loop and contribute to Acp significantly^. However, the experimental 



^ Lighter charged Higgs contribute dominantly to & ^ 57 and to satisfy that the constraint within the 3 a 
level one has to consider the non-minimal flavor violating MSSM scenarios. See Ref. 27[ for the details. 
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TABLE I: The maximal CP asymmetries and the number of required charged Higgs bosons for 
the four promising channels, H~ XiX2^ X2 XiJ X2 X-j) ^-^id H~ X2 X4- We also 

include the input parameter set for each channel. We show the results by including all experimental 
constraints and compare with the results where only mass bounds are consired. 



Chan. 


Acp{%) 


{Alp X Br)-\^) 


rriH- (GeV) 


MsusY(GeV) 


M2(GeV) 


/i(GeV) 


tan (3 


Consts. 


12 


-36.7 


3 X 10^ 


1040 


-600 


680 


-500 


41 


Mass Only 


12 


7.1 


1 X 10^ 


1320 


1040 


840 


-100 


44 


All 


21 


-24.5 


1 X 10^ 


1280 


640 


560 


460 


29 


Mass Only 


21 


-2.3 


1 X 10^ 


1440 


1000 


200 


220 


38 


All 


23 


-17.5 


3 X 10^ 


1120 


600 


240 


-460 


41 


Mass Only 


23 


-1.2 


6 X 10^ 


1160 


960 


200 


-500 


2 


All 


24 


-57.6 


2 X 10^ 


1480 


760 


720 


-180 


14 


Mass Only 


24 


8.8 


5 X 10^° 


1500 


960 


680 


-220 


5 


All 



constraints, especially EDMs, require a heavier SUSY scale, around 1 TeV. Note that we set 
the CP phase arg^Af) = 7r/2 for At of 1.2 TeV magnitude. For these values, mass parameters 
in the top scalar quark system affect the light quark EDMs (through chargino loops) and 
lighter SUSY scales are excluded to satisfy the EDM bounds. 

The most promising channel for observing a large asymmetry is XiX2^ which 

could show a maximal CP asymmetry around 7% with approximately 10^ x s^/e required 
charged Higgs produced at colliders. Even though H~ X2XA can produce a slightly larger 
asymmetry, it is not a better channel since it not only requires more Higgs particles, but also 
it is a decay channel into much heavier states. The maximal asymmetry for H~ X2X4 
requires 10^^ x s^/e events, much more than for the other three modes. Moreover, its 
branching ratio is of the order of 10~^. Thus the CP asymmetry in X1X2 has a 

better chance to be observed than the other three. Its branching ratio is around 2 x 10~^ 
which is small, but it is five orders of magnitude bigger than that of H~ X-zXl- Note 
that the distribution of the points with respect to zero axis is not symmetric and the CP 
asymmetry is mostly positive especially for the green points which satisfy all constraints 
(except H~ X2 Xa)- Though the relative sign is a prediction of MSSM, the absolute sign 
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is a result of choosing the phase of At to be positive. 

V. SUMMARY AND CONCLUSION 

We studied the CP asymmetry of the charged Higgs decays into pairs of charginos and 
neutrahnos in the framework of the MSSM. We neglected the possibility of flavor violation 
in the squark sector. We first analyzed the partial decay widths and the corresponding 
branching ratios for these channels. For a region of the parameter space with intermediate 
values of tan/5 and for ^ 600 GeV some chargino-neutralino channels are competitive 
with, and sometimes dominant over, other open channels. We then investigated the size of 
the asymmetry for various final decay states and the likelihood of observing it at colliders. 
A non-zero CP asymmetry requires an absorptive phase, for which we considered possible 
interference terms between tree-level and one-loop diagrams. We calculated the imaginary 
parts and presented both analytical expressions for the discontinuity in the diagrams with 
the help of Cutkosky rules, as well as numerical results. We have shown that the dom- 
inant phase is the one describing the left-right mixing the top scalar quark mass matrix 
and the asymmetry is maximum for aig{At) = 7r/2, while a.Tg{Ab) and aig{Ar) give negli- 
gible contributions. We analyzed the dependence of the asymmetry on the parameters of 
MSSM: MsusY, A^, tan/3, M2, [Af ;,,-|, arg(y4( ;,,-) and mn-- In each numerical and graphical 
investigation we include bounds from b — S7, {g — 2)^, Ap and the EDMs. 

We then performed a scan of the significant parameter space for maximizing the asym- 
metry, and correlate it to the number of charged Higgs needed to see the corresponding 
asymmetry shown to be proportional to {A'^p x Br)~^. In the scan we included the exper- 
imental bounds on the masses of supersymmetric particles only, as well as the previously 
mentioned low energy bounds. We first note that the inclusion of mass bounds restricts the 
asymmetry slightly, while the other bounds pose more severe restrictions on the CP asym- 
metries. In particular, the EDMs, thought to restrict parameters in the up and down scalar 
quark only, impose severe restrictions even on the masses and mixing of the third family of 
scalar quarks, confirming the expectation that the scalar quark mass parameters must be 
in the TeV region. Heavy scalar quark masses and heavy gluino masses are also needed to 
satisfy the 6 — * 57 constraint, while the bounds on — and Ap are not as stringent. Our 
analysis underscores the importance of including these constraints on the parameter space. 
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like in the analysis of top quark production and decays [27j. 

For the allowed regions, the most promising channel appears to be xf X25 "where 

both the asymmetry can reach 7% and the number of charged Higgs required to observe 
them is yet the smallest, of order 10^ x s^/e. This channel may be competitive with the 
quark channels Of course, this is encouraging, as light charginos and neutralinos are 
most likely to be produced in charged Higgs decays in the first place. By comparison, 
the other channels are less promising. For instance, for heavy chargino-neutralino channel, 
H~ — s> X2X^y almost 10^^ x s^/e charged Higgs events (with a branching ratio of the order 
of 10~^) are required for detecting an asymmetry of (8 - 9)%. 
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APPENDIX: THE METHOD - CALCULATION OF THE IMAGINARY PART 

We calculate the absorptive part of the loop diagrams by applying the Cutkosky rules. 
In general, in the loop integration, we end up with a numerator with scalar, vectorial, or 
tensorial structures (and their pseudo counterparts), depending on the types of particles 
running in the loop. For a more complete discussion of our method and the cut results for 



the two-point (self-energy) case and three-point (vertex-type) scalar case, please see [3|. We 
give here only the final formulas for the three-point (vertex-type) vectorial and tensorial 
cases with the horizontal cut: 



Ac>^ = ; ; ' log p^- 



{t - a^v) log 



Q3 + 
"3 - P3 



2/33 2pl^ , 



«3-/?3 V A 7 J 



VKpIpIpD 



^^^^^^2me{pl-{m, + m,y) \g^- _M..Ar^ APiPs+PiP'' 



zp^ zp^ 



2pl ^^PiP2+ ^^P3 
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AC 



1 



AG 



1 



as + P3 

«3 - P3 



2a3/53 



AC, 



1 



{3al - Pi + 4Aa3) log ( ) - 2f5^{?>a^ + 4A) 



«3 - 13^ 



+ 2/33 [as + Ai;(4t - 3a3f )] 



1(1 - 3At;2) + 2At(t - 2a^v)] log 



AC 



'P1P3 



AC 



2KpIpIpD 

+ 2/33[2t-i;(3a3 + 2A)]}, 

"3 + /?3 



[t;(3a^ - + 2a3(Aw -t)- 2\t\ log 



«3 + A 
"3 - A 



(t - U'iV) log 



2/53 



as -^3/ "3-/^3 



2r)2 

t-/53^ + ^(tt + 4p^) 



where 



«3_+A 
"3 - A 



. 2,2 2 

V = {pI-pI+pD/KpIpIpD, 

p3 = \/A(p?,pi,pi)A(pi,mf,mi), 

a = pUpI + 2m\ - (p^ + ^2 ^ + m^)] + (m^ - vr?^[p\ - pi), 
"3 = pH^pI + + P3 - + ml + 2ml)] + ("^1 - ^DiPi - pI)^ 
A = \{pi,pl,pl) 



(A.l) 



and (y'^'^ = (+, — , — , — ) is the metric tensor. Here mi and m2 represent the masses of the 
internal particles where we cut through. The four momenta p^* and P2/P3 momentum 
of the charged Higgs and neutralino/chargino, respectively. The other type of horizontal cut 
can be obtained by simply replacing ^3 ^ P2, and rrii ^ m2. 

For completeness we include the terms with the vertical cut, evaluated in 



AC: 



log 



27rz Q{pI — (mi + m2)^ 

VKpIpIpI) 

27ii 6(pi — (mi + m2)^) / g 



a + (3 
a — P 



Ps 



- ' (n - av) log ( + 2f3v^ 



a — P 



2 Pi 



pWKpl^phpl) 



AC, + ^ AC22 + (p^ Pl + p^) AC23 + P'3 Ps AC'ss 
z Hp^ , 
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Aa = —7^ - (a + /9) 
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